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In nuclear collisions at the LHC large number of hard partons are cre¬ 
ated in initial partonic interactions, so that it is reasonable to suppose 
that they do not thermalise immediately but deposit their energy and mo¬ 
mentum later into the evolving hot quark-gluon fluid. We show that this 
mechanism leads to contribution to flow anisotropies at all orders which are 
non-negligible and should be taken into account in realistic simulations. 


PACS numbers: 25.75.-q, 25.75.Ld 

Hot and dense strongly interacting matter which is produced in ultrarel- 
ativistic nuclear collisions like those at the LHC or RHIC, expands swiftly 
and soon decouples into individual hadrons. Ultimately, one wants to in¬ 
vestigate its evolution since it potentially bears an imprint of the Equation 
of State (EoS) and transport coefficients, which one wants to infer. 

When the distribution of produced hadrons in transverse plane is stud¬ 
ied, it appears far from isotropic even in most central collisions. Customar¬ 
ily, one studies its Eourier decomposition in azimuthal angle 
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where Vn and i/’n are parameters. Presently available statistics allows to 
determine the first six terms of such series even for individual events. 

The anisotropy is caused by anisotropic expansion of the hot fireball. 
Due to transverse expansion, transverse momentum spectra of hadrons are 
blue-shifted. This is because particles are produced by regions of the fire¬ 
ball which move transversely outwards in the direction towards the detector. 
Waves emitted by a source moving towards the detector are recorded with 
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shorter wavelength than emitted in the rest frame of the source. When 
translated in terms of momenta, this leads to an enhancement of higher 
momenta. This mechanism connects the transverse expansion of the fire¬ 
ball with the shape of the transverse momentum spectrum and maps the 
anisotropies of the former onto the anisotropies of the latter. The largest 
among the anisotropies is the so-called elliptic flow. It is mainly caused by 
the second-order anisotropy in the expansion velocity in non-central colli¬ 
sions due to initial anisotropic geometry of the fireball. 

A classical argument then connects direction dependent transverse ex¬ 
pansion velocities with inhomogeneities in pressure distributions within the 
fireball. Distribution of the energy density is established after the incoming 
partons interact and shows quantum fluctuations. This is in addition to 
any geometrical anisotropies which are due to non-zero impact parameter 
in the collision. The evolution from such an initial state is reasonably well 
described by hydrodynamic models. This is the point where EoS and trans¬ 
port coefficients enter the game. The resulting state of the expanded and 
cooled fireball depends on EoS and transport properties. Thus by calculat¬ 
ing hadronic spectra and their anisotropy one hopes to be able to tune e.g. 
the viscosities until an agreement with data is reached. 

There is a caveat, however, in this game. The initial conditions are 
not measured and can only be determined in various models. Unfortu¬ 
nately, they influence strongly the resulting anisotropies of the hadron dis¬ 
tributions [I]. Thus without the knowledge of the initial conditions the in¬ 
tended strategy for the extraction of transport properties seems jeopardised. 
Lucky enough, simulations with both ideal and viscous hydrodynamic mod¬ 
els show that there is quite a linear relation between the initial state spatial 
anisotropy and the anisotropy of hadron distribution Enam. This is true 
for the second and third orders and breaks for higher orders [ 2 ]. That can 
be understood, since higher order terms are influenced by the interference 
of lower orders (so that e.g. U4 gets contribution from the square of 62— 
second-order spatial distribution anisotropy). The linear relation allows to 
map the event-by-event fluctuations of V2 and U3 onto the fluctuations of £2 
and £3 and thus identify the model for the initial state which best agrees 
with the data. 

Usually, it is assumed that there is no contribution to fluctuations during 
the hydrodynamic evolution. This may not be the case, however. Hydro- 
dynamic simulation with the fluid energy and momentum density coupled 
to dynamically fluctuating order parameter field shows fluctuations in the 
energy density which can well cause flow anisotropies observable in data 
[ 5 ]. More precise quantitative impact of such a mechanism remains to be 
studied. 

Here, we introduce another mechanism that can induce flow anisotropies. 
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Hard partons from initial scatterings do not thermalise immediately as they 
are produced, but fly into the hot and dense quark-gluon plasma. There, 
mostly they are fully quenched so that their momentum is transferred into 
the fluid and must show up in the flow pattern. Directions of these partons 
are distributed isotropically, but due to their finite number there may be a 
contribution to flow anisotropy. 

In addition to that, the energy in a collision at the LHC is so large 
that there may be pairs of hard partons which are close in rapidity and are 
directed so that they might come close to each other during the evolution 
of the collision. Even if they are fully quenched before they could actually 
meet, their momentum has been shown to be further carried by the gen¬ 
erated streams in the fluid [H]. There is a good chance that such streams 
merge into one. In peripheral collisions, through this mechanism the collec¬ 
tive flow induced by hard partons tends to be directed in the reaction plane 
(which is spanned by the beam and the direction of the impact parameter). 

Imagine (Fig. the almond-shaped cross-sectional area of a fireball in 



Fig. 1. Transverse cross-section of the fireball in non-central collisions, with two 
dijet pairs produced. Reaction plane is horizontal. 

non-central collisions with two pairs of hard partons produced. If all four 
partons are directed parallel to the reaction plane, through the momentum 
loss they will all positively contribute to the elliptic flow (Fig. [^left). If, 
however, they are directed perpendicularly to the reaction plane, there is 
a chance that two of the streams will meet, merge and cancel flow in that 
direction (Fig. [fright). The chance is bigger in this direction because here 
the fireball is narrower and the streams have less space to avoid each other. 

We have set up a three-dimensional ideal hydrodynamic model and 
tested this scenario miH]. Note that a 3 D simulation in this case is manda¬ 
tory, since the presence of hard partons breaks the boost-invariance (and 
azimuthal symmetry). In order to esimate the effect on flow anisotropies 
due to hard partons we formulated the model without any other fluctua- 


4 


tomasik'eqcd printed on July 3, 2015 


tions in the initial conditions. Thus our initial energy density profile was 
determined from an optical Glauber model and the distribution in trans¬ 
verse plane was determined from a combination of wounded-nucleon and 
binary-collision density. The energy density in the central cell at the initial 
time of To = 0.5 fm/c was set to 60 GeV/fm^. The initial profile in rapidity 
is flat over 10 units with Gaussian tails at the edges. 

To implement the effect of energy and momentum deposition from hard 
partons into plasma, terms which represent forces are added into the 
energy and momentum conservation equation 
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where is energy-momentum tensor. The force field is parametrised with 
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where the terms in the bracket stand for the rate of energy and momentum 
deposition and the Gaussian {a = 0.3 fm) distributes it around the trajec¬ 
tory of the parton until all its energy is used up. The actual energy loss 
scales with the entropy density as 

dE _ dE 

dx dx 

with So = 78.2/fm^ (corresponding to Eq = 20 GeV/fm^) and dE/dx\Q 
being a parameter that we tuned. The places of hard parton production 
are distributed according to binary collision density in the transverse plane 
and uniformly in rapidity. Their directions are azimuthally symmetric and 
the pt spectrum follows from [9| 
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with B = 14.7 mb/GeV^, po = 6 GeV, and n = 9.5. 

In Fig. we show results obtained for collisions at vanishing impact 
parameter. There is no anisotropy in case of no hard partons as we initi¬ 
ate the simulation with smooth azimuthally symmetric energy density pro¬ 
file. Momentum deposition during the evolution leads to measurable flow 
anisotropies. Intristingly, the effect of energy loss does not depend on the 
value of dE/dx\Q. Note that the total deposited momentum for both tested 
values is the same, only the rate of deposition is changed. Simulations with 
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Fig. 2. Coefficients V 2 and V 3 from central collisions. Data calculated with two 
different values of hard parton energy loss, one simulation with only energy and no 
momentum deposition (hot spots). Shown also reference simulation with smooth 
initial conditions. 


only hot spots added in the initial state, which contain the same amount of 
energy as hard partons, but no momentum, lead just to a half of the effect 
of parton energy loss. 

We should prove yet, that the effect is correlated with event geometry so 
that in non-central collisions we indeed obtain positive contribution to the 
elliptic flow as argued above. To this end, events corresponding to impact 
parameter 6 = 6 fm (falls into the centrality class 30-40%) were simulated 
both with and without the hard partons contribution. Main results are sum¬ 
marised in Fig. where we show V 2 and V 3 as functions of pt- In non-central 
collisions, elliptic flow is generated due to the anisotropic matter distribu¬ 
tion from which the fluid evolves. Hard parton contribution increases the 
observed anisotropy by about 50%. As for U 3 , it is not present in case of 
smooth initial conditions due to missing third-order spatial anisotropy is 
solely generated from the hard parton momentum deposition. 

These results indicate that the contribution from hard parton momen¬ 
tum loss is significant and should be taken into account in hydrodynamic 
simulations that aim at the extraction of transport properties of the hot 
and dense nuclear matter. 

An open question is why changing the rate of energy loss has no influence 
on the generated elliptic flow. We speculate that the reason may be that 
practically all of the energy and momentum is deposited very early in the 
densest period of fireball evolution. We want to investigate this question 
more closely in the future. 

A realistic simulation, also to be accomplished in the future, must in¬ 
clude viscosity effects and fluctuating initial conditions. There we plan to 
use a newly developed hydrodynamic model m- 
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Fig. 3. Coefficients V 2 and V 3 from collisions at impact parameter 6 = 6 fm. We 
show results from simulations with only smooth initial energy density profile (leg¬ 
end: no partons) as well as with energy loss from hard partons included. 
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